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An approximate local potential for the residual tt^-tt^ interaction is computed. We use an 0{a^) improved 
action on a coarse 9^ x 13 lattice with a « 0.4fm. The results present a continuation of previous work; Increasing 
the number of gauge configurations and quark propagators we attempt extrapolation of the tt^-tt"'" potential to 
the chiral limit. 
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1. INTRODUCTION 

The use of improved lattice actions allows to 
work with lattice volumes large enough to accom- 
modate systems of two hadrons, with manageable 
computational effort. We take advantage of this 
opportunity to study the residual effective inter- 
action of two pseudo-scalar mesons on the lattice. 

This paper is a report on the current status 
of this project. An earlier exploratory study [n| 
was made with a smaller number of gauge config- 
urations and two, somewhat large, quark masses. 
The new results presented here are based on 208 
configurations and 6 values of quark nrasses, with 
otherwise unchanged lattice parameters. The 
current simulation allows extrapolation of the 
extracted tt-tt potential to the chiral limit and 
a comparison of lattice-based scattering phase 
shifts, computed with the potential, to experi- 
mental results [0,^ in the isospin 1 = 2 channel. 

2. LATTICE PARAMETERS 

An L"^ X T = 9^^ X 13 lattice was used with 
an 0{a'^) tree- level and tadpole improved ac- 
tion based on next-nearest neighbor couplings 
|4|§. At /3 = 6.2, in the conventions of [|, 
the corresponding lattice constant is a w 0.4fm, 
or a-i « 500MeV. We have used Nu = 208 
quenched gauge configurations. The hopping pa- 



rameters for Wilson fermions were set to k ^ = 
5.720,5.804,5.888,5.972,6.056,6.140. The criti- 
cal value for k~^ is « 5.5. Quark propagator 
matrix elements were computed using a randonr- 
source technique with Nj^ = 8 Gaussian 
sources. 

3. METHOD 

An outline of the theoretical framework may be 
found in pp] , we here mention only the essential 
points. Suitable operators for the tt+-tt^ system, 
having isospin 1 = 2, are 

^p{t) = 0_j7(i) 0+p-(i) , (1) 

where 

Mt) = L'^ E ^'^'""^"{^^ t) i^rix, t) (2) 

X 

describes single mesons with lattice momenta p. 
The correlation matrix for the tt^-tt^ system 



Cp,- = {^li^ *g(io)) = Cp^+Ci,f;^ 



(3) 



is a sum of a free, C, and a residual-interaction 
contribution, Cj. The free TT+-7r+ correlator C is 
diagonal in p, q. We also implement link variable 
fuzzing and operator smearing |||,y at the sink. 
We define an effective interaction through 
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iJj = _^ln((7-l/2c C-l/2). 



(4) 



Matrix elements of Hi are obtained from linear 
fits to the logarithm of the eigenvalues of the cor- 
relators C and C At this point, only the diag- 
onal elements of C were utilized. Momentum- 
space matrix elements {p\Hi\q) are computed in 
a truncated basis of small lattice momenta. The 
Fourier transform to coordinate space contains a 
local potential 



v{f)^j2' 



— 2iq-r 



{-q\Hi\ + q: 



(5) 



Its s-wave projection {(. — 0) makes use of only 
the irreducible representation Ai of the lattice 
symmetry group 0{2>,Z). In terms of the cor- 
responding reduced matrix elements we have 

FW=^Jo(2<zr)(g|ff(^^)|<z), (6) 

where jo is a spherical Bessel function. 

4. RESULTS 

Potentials according to (p|) for one selected 
value of K, are shown in Fig. nT The sum over (on- 
axis) momenta q = ^k was truncated at increas- 
ing fcniax — 0,1,2,3,4, respectively. A detailed 
error analysis is in progress, however, errors ap- 
pear to increase with k. The results presented 
below are for the truncation at fcmax — 3. 

Chiral extrapolation of the potential was done 
by linear fits of V{r) versus m^ for a fine (plot- 
grade) mesh of values r, fixed one at a time. Us- 
ing sets of 3 through 6 data points, corresponding 
to the smallest available values of m^, gives very 
similar results. The subsequent analysis was per- 
formed with 3 data points. 

The extrapolated potential V{r) is shown in 
Fig. H as a dashed line. The oscillations of V{r) in 
the region r > 2 are due to the Fourier transform 
of the truncated momentum sums. The wave 
length is indicative of the lattice resolution at the 
current truncation. 

A parametric fit to V{r) with 
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at a5 = 2 fixed, was applied to the extrapolated 
potential. The result is shown in Fig. as a solid 



Figure 1. Local potentials (|^) for the residual tt-tt 
interaction at k"^ = 5.720. Results using trun- 
cated correlation matrices of sizes one-through- 
five (fcinax — 0, 1, 2, 3, 4) give a hint at the conver- 
gence behavior of the Fourier series. 



line. It suggests attraction at short distances fol- 
lowed by a repulsive barrier. 

We have used V'-'^-^ir) in a Schrodinger equa- 
tion to calculate s-wave scattering phase shifts 
SJ,^q{p), see Fig. ||. The pion mass was set to 
multiples of the experimental value, correspond- 
ing to rriT; = 0.28 in units of a^^. The repulsive 
nature of the phase shifts is due to the hump of 
V^'^'jr) around, and extending beyond, r w 1, see 
Fig. 0. The data points in Fig. are experimental 
results compiled from H,^. 

5. CONCLUSION 

Scattering phase shifts for the 1 — 2 channel 
TT-TT system were computed from lattice QCD by 
way of extracting a non-relativistic potential. 

Since at this point an error analysis is pend- 
ing, particularly systematic, errors are unknown. 
The range of the extracted potential is short com- 
pared to the current spatial resolution. The latter 
is determined by the somewhat large value of the 
lattice constant, and by the limitations imposed 
by the momentum truncation for the correlator 
matrices. This situation makes it difhcult to reli- 
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Figure 2. Result F(r') of the chiral extrapolation 
shown by the dashed line. A parametric fit with 
^("^(r), see (0), gives the solid line. 



ably extract details of V{r). (The current lattice 
parameters should be better suited for studying 
interactions involving, larger sized, baryons.) In 
addition to V{r) there is also present a nonlocal 
potential (l) which has not yet been computed. 

The scattering phase shifts obtained from the 
underlying lattice study show repulsive behav- 
ior in the low-momentum region and, in this re- 
spect, compare favorably to experimental find- 
ings. Quantitatively, the computed phase shifts 
are too small by a sizeable factor. Relativistic 
corrections are at the 40% level at p « 0.6a~^. 
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